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Abstract 

Two-photon annihilation rates of 2"*" tensor charmonium and bottomonium up to third radial 
excited states are estimated in the relativistic Salpeter method. Full Salpeter equation for 2+ tensor 
state is solved with a well defined relativistic wave function and we calculated the annihilation 
amplitude using the Mandelstam formalism. Our estimates of the decay widths are: T{xc2 
27) = 501 eV, r(xc2 ^ 27) = 534 eV, T{xb2 ^ 27) = 7.4 eV and T{x'b2 ^ 27) = 7.7 eV. We also 
give total decay widths of the lowest states estimated by the two-gluon decay rates, and the results 
are: Ttot{Xc2) = 2.64 MeV, Ttot{xb2) = 0.220 MeV and Ttot{Xb2) = 0-248 MeV. 
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I. INTRODUCTION 



Recently, the radiative annihilation physics of Xco, Xbo (0'*'^), Xc2 and Xb2 (2^^) become 
hot topics [l|, y, y, y, y, y, y, is, y, l^, hi, because the annihilation amplitudes are related 
to the behavior of the wave functions, so the annihilation rates are helpful to understand 
;he formalism of inter-quark interactions, and can be a sensitive test of the potential model 



121. 



In previous letter [13|], two-photon and two-gluon annihilation rates of 0"''+ scalar cc 
and bb states are computed in the relativistic Salpeter method, good agreement of our 
predictions with other theoretical calculations and available experimental data is found. In 
our calculation, we found the relativistic corrections are large and can not be ignored, and 
point out that all the calculations related to a P wave state, one have to use a relativistic 
method, a „on-.eMiv.t. method „m cau. a la.ge e™. even fo. a hea^ =tate fl Q. 
In a non-relativistic calculation, the corresponding decay width is related to the derivative 
of the non-relativistic P-wave function at the origin, but in a full relativistic calculation, 
the relativistic corrections include not only relativistic kinematics but also the relativistic 



inter-quark dynamics, the decay width is re 



which can be seen in this letter or in Ref. 



15| 



ated to the full behavior of P-wave function 



In this letter, we give relativistic calculation of 2+"*" tensor cc and bb states decaying into 
two photons usin g th e instantaneous Bethe-Salpeter method pi3j] , which is a full relativistic 
Salpeter method [17|. The case of tensor 2"*"+ state is special, not like other P-wave states, 
there are no decay constants in this state, and because there is the P — F mixing problem, 
it make the physics much complicated. 

The reason of the exist of P — F mixing is that P wave and F wave state have the 
same parity and charge conjugate parity, they are all J^*" = 2+"'', one can not distinguish 
;hem by the quantum number. We have found a basic method to deal with this problem 
isl . we begin from the quantum field theory, analyze the parity and charge conjugation of 
bound state, and give a formula for the wave function that is in a relativistic form with 
definite parity and charge conjugation symmetry, then we solve the full Salpeter equation, 
and obtain the corresponding state, and there are automatically the mixing between P wave 
and F wave. 

The letter is organize as following. In Sec. II, we give the annihilation amplitude in Man- 
delstam formalism and the wave function of the 2"*" tensor state with a well defined relativistic 
form. The two-photon decay width and full width of heavy 2^^ quarkonium are formulated 
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FIG. 1: Two-photon annihilation diagrams of the quarkonium. 



in Sec. II, we show the numerical results and give discussions in the Sec. III. 



II. THEORETICAL DETAILS 



According to the Mandelstam [19| formalism, the relativistic transition amplitude of a 
quarkonium decaying into two photons (see figure 1) can be written as: 



27 



{2ny 



hS{pi - h) ^1 + ^iS{pi - k2) ^2 



(1) 



where ki, ei, 62 are the momenta and polarization vectors of photon 1 and photon 
2; = I for charm quark and Cg = | for bottom quark; pi and p2 are the momentum 
of constitute quark 1 and antiquark 2; is the Bethe-Salpeter wave function of the 

corresponding meson with the total momentum P and relative momentum q, related by 



p^ = aiP + g, «! 
P2 = a2P - q, a2 



mi 



mi + 1712 

1712 



Till + m2 

where mi = m2 is the constitute quark mass of charm or bottom. 

Since pio + pp.o = M, the approximation piQ = p2o = ^ is a good choice for the equal 



mass system 



211 ]. Having this approximation, we can perform the integration over go 
to reduce the expression, with the notation of Salpeter wave function \l/(g) = / ^xl?)? to 



T2^ = Vs (ee^ 



1 7^. tr { m 



7 h + h -j- 

k - mi pi- 



(2) 



t2 - rrii 

The general form for the relativistic wave function of tensor = 2+ state (or J^^ = 2+"*" 
for equal mass system) can be written as 16 terms constructed by momentum P, q and Dirac 
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matrix 7, because of the approximation of instantaneous, 8 terms with P ■ q become zero, 
the relativistic Salpeter wave function for 2+ state can be written as: 



/i(g) + ^/2(g) + §/3(g1 



M2 



(3) 



where the Ef^i, is the polarization tensor of the 2"*" state, q± = (0, g). But these 8 terms wave 
unctions fi are not independent, there are the further constraint from Salpeter equation 



171]: "$2+ (q) — ^2+^(^ ~ ^ ' which give the constraints on the components of the wave 



function, so we get the relations 

[qlhiq) + M'Mq}] (uj, + UJ^) - M^/sl^K^i - ^^2) 



f2{q) 



M{miuj2 + m2UJi) 
[qlfM'M^U{q)]{u,~U2) 



M{miU2 + m2U!i) 



U{q)M{uj^-uj2) 



Uq)M{uJi + tU2) 



(4) 



miLJ2 + 1712001 miLJ2 + m2U!l 

Only four independent wave functions /3(g), /4(g), /5(g) and feiq) been left, one can check in 
Eq.(l3]), all the terms except the two terms with /2 and are charge conjugate parity positive, 
but /2 and terms have negative charge conjugate parity, after we use the constraint 
relations, for equal mass system, the terms with /2 and fj become zero, then the whole wave 
function have positive charge conjugate parity, that is 2++ state. These wave functions and 
the bound state mass M can be obtained by solving the full Salpeter equation with the 
constituent quark mass as input, and we will not show the details of how to solve it, we only 
show our result in next section. 

These four independent wave functions fulfil the normalization condition: 



dq 



16 UJ1UJ2 (f 



(27r)^ 15(miCo'2 + ^21^1 

(mi + m2){m2UJi - miU2 



5 + 



{nil + m2){m2UJi - miU2) 

LL)iLL)2{i^l + ^^2) 



7^ 



2 + 



2g73 /4^ + /6 



2M. 



(5) 



UJiUJ2{uJl + UJ2) 

With the full Salpeter wave function, the two photon decay amplitude can be written as: 

dq \ 1 



T, 



27 



4 ^3 £ 



(27r)3 \ xi-\q\ M cos9 



[ h M (A;f ql ■ 62 



+6^ ql £2 ■ g± + 4 ql ^1 ■ qi-) + ^ {ei ■ £2 g± ■ + 2 ei ■ £2 ■ g±) 
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Xi + \q\ M cosO 



[ h M {kl^ ql £i ■ £2 + ql £2 ■ g± + 4 ql ^1 ■ g±) 



+ ^ ^ {£1 ■ 62 q±- k2 + 2 El - q^ 62- q±) 



(6) 



where Xi = ^ + q + ml, and 9 is the angle between the momentum q and ki. Finally 
the decay width with first order QCD correction ^] can be written as: 

= 2!.5-16 ttM ^ ■ " ^) • 

Until now, only the total decay width of Xc2{^P) is available, we can estimate the full 
decay width of OZI-forbidden states using the two gluon decay width, and the two gluon 
decay width of quarkonium can be easily obtained from the two photon decay width with 
a simple replacement. For the charmnium states, only the ground state Xc2(l-P) is OZI- 
forbidden, we have: 

^ 9 a^(mc) 1 - ""^""'^^ 

^totiXc2) = r2g(Xc2) = r2^(Xc2) g' o ^ . 16 Mm,)" " i^) 

3 TT 

For 2~^~^ bottomnium states, according to our estimate of mass spectra, there are two 
states which are below the threshold of BB, and we can predict their full decay widths 
using their two gluon decay widths. For Xb2{^P), we have: 

rtot{xb2) = r2,{xb2) = r2M Zjm,) > (9) 

and for the first radial excited state xJ,2(2-P), we have: 

18 a'^{mi)) 1 + 

^tot{Xh2) - ^2g{Xh2) = r27(Xfe2) 16 ^m,)" ' i^^) 

3 TT 



where a = and the QCD corrections are summarized in Ref.[22|. 



III. NUMERICAL RESULTS AND DISCUSSIONS 

We will not show the details of Solving the full Salpeter equation, only give the final 



results, interested readers can find the the detail technique in Ref. 23 1 



When solving the full Salpeter equation, we choose a phenomenological Cornell potential, 

ng)=K(g1+7o®7°K(g-), 
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Vsiq) = -i- + Vo)5'{q) ^ 



37r2 (g^ + a2 

and the coupling constant as{q) is running: 

127r 1 



^'^(^1 = -£i7^^ ' (11) 



25 log(a + ^ 



^QCD' 



Here the constants A, a, a, Vq and Aqcd are the parameters that characterize the potential 
We found the following best-fit values of input parameters which were obtained by fitting 
the mass spectra for 2++ Xc2- a = e = 2.7183, a = 0.06 GeV, Vo = -0.401 GeV, A = 0.2 
GeV^, Aqcd = 0.26 GeV and rric = 1.7553 GeV. With this parameter set, we solve the full 
Salpeter equation and obtain the mass spectra shown in Table I. To give the numerical value, 
we need to fix the value of the renormalization scale /z in a<i(/i). In the case of charmonium, 
we choose the charm quark mass rric as the energy scale and obtain the coupling constant 



a,(mj = 0.36 |23|. 



Since the wave function Eq.([3]) is general wave function for state, and either and 
■^^2 can be 2+ state, so the obtained states of x'c2i^^) Xc2(2-^) are not pure F wave, 
they are P — F mixing state, but a F wave dominate state, and Xbii^^)-! x'LC^^) Table 
II are also P — F mixing, but F wave dominate state, we will discuss the detail of mixing 
in other paper 24 1 . 

Our prediction of the mass for x'^2(2-P) is 3967.0 MeV, which is a little larger than the 
first observation by BELLE collaboration, their data is 3931 MeV. And our prediction of the 
first F wave dominate state, we list as x'c2i^F)^ whose mass is 4040.5 MeV, and the second 
one, Xc2(2-^), whose mass is 4314.3 MeV. 

With the obtained wave function and Eq.Q, we calculate the two-photon decay width 
of cc 2++ states, the results are also shown in Table I. Not similar to the S wave case, where 
the two photon decay width of ground state T{ric 27) is much larger then the first radial 
excited state T{rj'^ 27), almost twice p^], the ground state decay width is a little smaller 
than the first radial excited state, and from the table we obtain the conclusion that the 
decay widths with successive radial excitation fall very slowly. The decay width of the F 
wave dominate state, T{x^2{^P) ^ 27) = 92 eV and r(Xc2(2-F) 27) = 54 eV are much 
smaller than the P wave dominate state. 

For the case of 66, our best fitting parameters are Vq = —0.459 GeV, Aqcd = 0.20 and 
nib = 5.13 GeV, other parameters are same as in the case of cc. With this set of parameters. 
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the coupling constant at scale of bottom quark mass is asirrib) = 0.232. The corresponding 
mass spectra, two-photon decay widths are shown in Table II. Our mass prediction of x'b2{'^-P) 
is 10283 MeV, a little higher than the data. And the two photon decay widths are much 
smaller than the case of charmonium, their value is only about several eV. 



TABLE I: Two-photon decay width and total width of P-wave 2'^~^ charmonium states, where the 
data of Xc2(l-P) is come from PDg[25'], and data of Xc2(2-P) is come from Ref.[26']. 





Ex Mass (MeV) 


Th Mass (MeV) 


T2^ (eV) 


Xc2{lP) 


3556.20 


3556.6 


501 


Xc2m 


3931 


3967.0 


534 


x'c2m 




4040.5 


92.4 


Xc2m 




4264.6 


509 


x:2m 




4314.3 


54.1 


xlm 




4498.7 


475 


TABLE IL Two-photon decay width and total width of P-wave 2^^^ bottomonium states. 




Ex Mass (MeV) 


Th Mass (MeV) 


T2^ (eV) 


Xb2{lP) 


9912.21 


9912.2 


7.43 


xmm 


10268.65 


10283 


7.69 


Xb2m 




10364 


1.76 


xl2m 




10561 


7.19 


Xmim 




10616 


1.43 


x;;(4P) 




10786 


6.60 


Our predictions of the total decay 


width for the ground state: 






^tot{Xc2) 


= r29(Xc2) = 2.64 MeV, 




is little larg 


;er then the PDG data Ttot 


(;^^2) = 2.03 ±0.12 MeV. For 2+ 


+ bottomnium states. 


we have no 


data in hand, our theoretical predictions are: 






^tot{Xb2) 


= r2,(Xfe2) = 0.220 MeV 




and 










'^tot{x'b2) 


= r29(X62) = 0.248 MeV. 
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TABLE III: Recent theoretical and experimental results of two-photon decay width and total width. 



rgMkcV^ r^„f (MeV) rl;^ (kcV) r^^^ (eV) r^^^ (MeV) r^^^ (eV) 



Ours 0.50 2.64 0.53 7.4 0.22 7.7 

Gupta [6j 0.57 1.20 8 0.22 

Ebert [8] 0.50 0.52 8 6 

Miinz [9] 0.44±0.14 0.48±0.16 5.6±0.6 6.8±1.0 

CLEO [27| 0.53±0.15±0.06±0.22 

CLEO [28| 0.56±0.06±0.05±0.04 

CLEO [29] 0.60±0.06±0.06 

PDG [251 0.493 2.03±0.12 



We compare our predictions with recent other theoretical relativistic calculations and 
experimental results in Table III. Except the total decay width of Ttot{Xc2)i all the values of 
listed in the table consist with each other. 

In summary, by solving the relativistic full Salpeter equation with a well defined 
form of wave function, we estimate two-photon decay rates: T{xc2 27) = 501 
eV, r(x;2 ^ 27) = 534 eV, r(xfe2 27) = 7.4 eV and r(x62 ^ 27) = 7.7 eV, and the 
total decay widths: Ttot{Xc2) = 2.64 MeV, Ttot{Xb2) = 0.220 MeV and Ttot{Xb2) = 0-248 MeV. 
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